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Silicon Burning and the Synthesis of Cosmic-Ray Nuclei 

Donald V. Reames 

NASA/Goddard Space Flight Center 
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ABSTRACT 

The source abundances of the cosmic-ray nuclei of atomic 
number 14 < Z • 28 are found to be in good agreement with 
calculations based on the quasi-equilibrium silicon-burning 
process. The low-density environment appropriate for cosmic- 
ray nucleosynthesis differs from that required for synthesis 
of the observed solar-system ("universal") abundances in a way 
which is consistent with the shock acceleration model for the 
supernova origin of the nuclei observed in both abundance 


measurements 


Tho quasi-equilibrium silicon-burning mechanism has been success- 

1,2 

fully applied to the explanation of the solar abundances of the 
nuclei with atomic numbers in the range 14 • 2 ' 28. This mechanism 
operates when Si ' in a free bath of alpha particles is heated to 
temperatures above 3 x 10 9 K, with the heavier element abundances 
being linked to that of Si r ' primarily through a chain of (a, v ) reactions 
that are in equilibrium with the inverse reactions. At 4,4 x 10 y K and 
silicon is rapidly burned with a time scale of «' 1 sec; these time 
scales become compatible with that of tho shock wave in a supernova^ 
where the mechanism is believed to operate^. 

Supernovae have been recognized as possible sources of cosmic 

5 3 

radiation for some time and the model studied by Colgate and White 

describes the acceleration of the outer layers of the supernova to 

relativistic energies. It, therefore, seems reasonable to attempt to 

ascribe the composition of cosmic radiation in the region from silicon 

through iron to silicon burning e The cosmic-ray composition is known 

to differ from that of the sun^ even after corrections have been made 

for the passage of the radiation through interstellar material. In the 

following we find that those differences are not inconsistent with 

silicon-burning in an environment of lower density (o ~ lCf g/cm ,? ) than 

that required to obtain agreement with solar abundance (o * KTg/cm 3 ). 

This difference in density is completely consistent with the supernova 

3 

shock-acceleration model . 


2 



If we confine our attention to cosmic-ray nuclei at relativistic 

energies where changes in energy are small, the flux J^(x) of i-type 

nuclei above any given energy after passing through x g/cnf of inter- 

5 7 

stellar material is described by * 

dj. (x)/dx *-Ji(x)/\j+ J k (x)A ik (1) 

1 k>i 

where X^ is the mean free path for loss of i-type nuclei and X^ k is the 
mean free path for production of i-type nuclei from (heavier) k-type 
nuclei. 

5 

It has been found that the distribution of matter traversed is 

well represented by an exponential so that the observed flux is 

found (within a constant factor) from 

. 00 

jf ■ ■ Ji(x) exp (-x/x 0 ) dx ( 2 ) 

This distribution is appropriate when a time equilibrium is established 
between the loss of cosmic rays by galactic escape (with mean free path 
x 0 ) and their production from a uniform continuous distribution of 
sources, but it is also approximately valid for a random distribution 
of discrete sources distributed within the galactic disc® 9 ^. 

If we note that is formally the Laplace transform with respect 
to l/x Q of J^(x) and take the Laplace transform of Eq. (1) we find 


so that the flux of i-type nuclei from the source is a linear combi- 
nation of observed fluxes* Eq.(3)may also be obtained as the energy 
independent limit of the solutions of the equilibrium problem^***. 

We have applied Eq. (3) to the observed cosmic-ray composition 
summarized by Shapiro and Silberberg^. The fragmentation mean-free- 
paths were calculated using the empirical cross-section formula of 
Rudstam^ and that of Shapiro and Silberberg^. Differences between 
the results using the two formulae were not found to be significant 
for the nuclei of interest and the Rudstam formula was used in obtain- 
ing the cosmic-ray source composition shown by circles in Fig* 1. The 
observed composition used for these calculations is shown by crosses 
in the figure. 

The value of x Q is customarily determined by the requirement that 

Li, Be and B are purely secondary nuclei. Using Eq. (3) and the frag- 

13 

mentation cross-sections of Yiou , this condition occurs for 

~ ^ gm/cm 2 . For this value, however, the source fluxes Cl, K, Sc, 

Ti and V become negative, implying that the observed fluxes should be 

increased by as much as three standard deviations for K and Sc. If 

we require that these fluxes be increased by no more than two standard 

2 

deviations we find x Q £ 4 g/cm in calculating the results shown in 

Fig* 1. This lower value of x is also in better agreement with the 

o 

composition of the nuclei with 30 ^ Z ^ 90.^ 

The composition resulting from Si burning has been calculated 

p 

from the results tabulated by Bodansky, Clayton and Fowler*. Within 
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the range of environments studied by those authors, 3,6 x 10 J T 
• 5.0 x 10 <J K and ltf’ * p ■. 10 y g/cm M best agreement with the cosmic- 
ray data is obtained for T * 4,6 x 10 u K, p * irf'g/cm*' with the 
ftaction of silicon remaining, f * .15; these results are shown in 
Fig. 1. The time scale for the burning is t * 77 msec, for this 
case. 

For any temperature and density the comparable abundances of Fe 
and Si imply 0 15 ^ f < .25, However, the low abundances of elements 
in the region 15 • Z £ 19, a distinctive feature of the cosmic 
radiation, are obtained only for the lowest densities at a given 
temperature. The results are not highly sensitive to temperature 
within the above constraints. Since the time scale is a strong 
function of temperature, it is also poorly determined. 

Calculations with p > 10 P g/cnf* , the value used for comparison 
with the solar composition^, result in abundances of S, Ar and Ca which 

are factors of 2 , 3 and 6 larger, respectively, than those for 
p = 10 L g/cm 3 shown in the figure, liven at 10 5 g/cnf the calculated 
sulphur abundance is more than one standard deviation above the 
cosmic-ray abundance suggesting that even lower densities might be 
appropriate. 

For a supernova of about two solar masses, the theory of Colgate 

3 

and White predicts that an external mass fraction of ~10" 4 will be 
ejected as relativistic cosmic rays. The environment attained in these 
outer shells of the supernova is in excellent agreement with that which 
we find from the cosmic-ray composition. 
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We note that material burned rapidly at low density seems in- 
appropriate to a post- supernova environment which would obtain in 
the vicinity of a pulsar. 

A significant feature of the low density solution is that Ni bii 
can dissociate into Fe b4 + 2p in this environment which is neutron 
rich relative to the solutions at higher density, free neutron to 
proton ratio being 1.7 x 10" 4 herv. ^vs 0 ~ 10" 8 at p * 10 H g/cm 3 ). 

Thus the dominant iron isotope produced is Fe 34 rather than Fe b(X 
resulting from the decay of Ni DtJ which is produced deeper in the super- 
nova envelope. Since we have not attempted to exclude lower density 
solutions which might result in the direct production of Fe Bs , the isotope 

structure of cosmic-ray Fe is not clearly determined. A large pro- 

F 15 

portion of Fe 4 in cosmic rays would affect the recently-suggested 

technique for determining the age of the radiation by observing the 

effects of Mn 53 , though not adversely since Mn 63 would be the only 

quasi-stable isotope of Mn produced as a secondary from Fe 64 fragmenta* 

tion in interstellar hydrogen. 

Insofar as the nuclei in the range 14 ri Z < 28 are concerned it 
seems possible to form a consistent model for their synthesis via the 
quasi-equilibrium silicon-burning mechanism in supernovae. The bulk 
of the material ejected is burned at high density and later accreted 
by stars such as the sun while the outer layers, burned at lower density, 
are accelerated to high energies and are observed as cosmic radiation. 

1 

Postulating that the same mechanism is responsible for both samples of 
material clearly does not imply the same composition for both samples. 
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In fact the model would imply gradual changes in the cosmic- ray 
composition with energy since cosmic rays accelerated to different 
energies would experience nucleosynthesis in different environments. 

The observation of changes in cosmic-ray composition with energy 
eventually could provide a powerful experimental test of the 
mode 1 . 

It should clearly be emphasized that the present work is con- 
fined to a very limited region of nuclear charge. The abundances of 
other nuclei in cosmic radiation, especially those with 6 < Z •* 14, 
could be highly relevant and remain unexplained by the present model. 
However, this being the first known attempt at a quantitative comparison 
of nucleosynthesis theory with cosmic-ray observations, v/e find the 
agreement obtained here encouraging. 
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FIGURE CAPTION 


The abundance of elements extrapolated to the cosmic ray 
sources (solid circles) arc compared with the results of 
silicon burning with T * 4.6 x 30° K, p * lO^g/cm 5 and 

f ■ ,15 (see text). Also shown are the observed cosmic-ray 
abundances (crosses). All abundances are normalized to 


that of silicon. 
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